Phases of the infinite U Hubbard model 
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We apply the density matrix renormalization group (DMRG) to study the phase diagram of the 
infinite U Hubbard model on 2-, 4-, and 6-leg ladders. Where the results are largely insensitive 
to the ladder width, we consider the results representative of the 2D square lattice model. We 
find a fully polarized ferromagnetic Fermi liquid phase when n, the density of electrons per site, 
is in the range 1 > n > uf ~ 4/5. For n = 3/4 we find an unexpected commensurate insulating 
"checkerboard" phase with coexisting bond density order with 4 sites per unit cell and block spin 
antiferromagnetic order with 8 sites per unit cell. For 3/4 > n, the wider ladders have unpolarized 
groundstates, which is suggestive that the same is true in 2D. 



-t— > 

CO 



I 

C 

O 

o 



> 
in 

cn 
rn 

o 



The Hubbard model is the paradigmatic representa- 
tion of strongly correlated electron systems in condensed 
matter physics. It is widely invoked in studies of metallic 
ferromagnetism, unconventional superconductivity, vari- 
ous forms of charge and spin density wave formation, and 
even in theoretical explorations of spin-liquid and non- 
Fermi liquid states. These studies often involve conflict- 
ing claims concerning the phase diagram of this model 
in more than one dimension, especially when the inter- 
action strength U is comparable to or larger than the 
bandwidth. 

In this paper we report the results of an extensive 
DMRG study of the zero temperature (T = 0) phase 
diagram of the Hubbard model in the U — > oo limit, as a 
function of n, the density of electrons per site. To begin 
with, we study the two leg ladder on systems of size up 
to 2 x 50, large enough that finite size scaling can be 
used to obtain clear convergence to the thermodynamic 
limit. The resulting phase diagram is shown in the lower 
panel of Fig. 1. It is frequently found that important 
features of two dimensional (2D) quantum models are al- 
ready apparent in the solution of the 2-leg ladder. To get 
a feeling for which features of the 2-leg Hubbard ladder 
extrapolate smoothly to 2D, we compute the properties 
of 4-leg and 6-leg ladders (with sizes up to 4 x 20 and 6 
x 8, respectively). The inferred partial phase diagrams 
of these wider ladders are shown in the two middle panels 
of Fig. 1. (We have also carried out limited additional 
studies of 3 and 5 leg ladders.) While there may be sub- 
tle correlations characteristic of the 2D model that would 
only be manifest were we able to study wider or longer 
ladders, many features of the phase diagram are already 
remarkably insensitive to ladder width and length for the 
studied system sizes. We therefore speculate that these 
features survive as groundstate phases of the fully 2D 
model, as shown in the upper panel of Fig. 1. 

To summarize our findings: 

1) For 1 > n > rip, we find a half-metallic ferromag- 
netic (HMF) phase, i.e. a fully spin polarized Fermi liq- 



uid. For all even leg ladders we have studied uf = 4/5, so 
we expect that in 2D, rip ~ 4/5 as well. As is well known, 
the existence of a HMF phase for n close enough to 1 is 
suggested by Nagaoka's theorem, [2] and by various previ- 
ous exact diagonalization[3] and variational [4, 5] studies, 
the most recent of which [5] yields estimates of np which 
are more or less consistent with the present results. On 
the other hand, other lines of analysis [1] are suggestive 
that the HMF is stable only below a critical "doped hole 
density," (1 — n), which vanishes in the thermodynamic 
limit (np — ¥ 1~); the present results falsify this. 

2) For uf > n > 3/4, the 2-leg ladder appears to phase 
separate (PS), with the two coexisting phases having den- 
sities n = tif and n — 3/4. Our more limited results on 
broader ladders suggest that the same might hold true 
for 4 leg ladders and, by extension, in 2D as well. 

3) For n = 3/4, the 2-leg ladder forms an insulating 
commensurate plaquette density wave state, as shown 
in Fig. 2 (a). This pattern of symmetry breaking sug- 
gests that the spin-degrees of freedom correspond to an 
effective spin 3/2 Heisenberg antiferromagnetic chain on 
the "checkerbord" lattice defined by the "strong" plaque- 
ttes; indeed there is no detectable spin gap and we find 
clear signatures of quasi-long-range antiferromagnetic or- 
der with twice the period of the plaquette order, as also 
illustrated in Fig. 2(a). The 4-leg ladder exhibits a sim- 
ilar (and slighly weaker) ordering tendency, forming the 
checkerboard plaquette order of the sort shown in Fig. 
2(b); the spin correlations are consistent with those of an 
effective 2-leg spin 3/2 antiferromagnet on the checker- 
board lattice in that there appears to be a small spin-gap, 
with the implication that, in this case, the correspond- 
ing antiferromagnetic order (also shown in the figure) 
decays exponentially at long distances. Given that the 
6-leg ladder (Fig. 2(c)) also shows a similar ordering pat- 
tern as the 4 and 2-leg ladders, we arc confident that the 
corresponding phase persists in the 2D limit. The 2D 
"checkerboard" phase has coexisting bond-density wave 
and block-spin antiferromagnetic order [9]. However, the 
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site-charge density is uniform. The existence of this 
phase, and its apparent robustness, was unanticipated 
in previous studies, as far as we know. 

4) For 3/4 > n > n p with n p w 3/5, the 2-leg ladder 
has a partially polarized ferromagnetic ground-state, 
with maximal spin-polarization for n ks 2/3. We have not 
yet determined where, in this interval, the system forms 
a phase separated two-phase mixture and where it is in a 
pure phase. In any case, whatever the phase diagram of 
the 2-leg ladder, it is unlikely to be representative of the 
2D system in this range of n, given that the groundstate 
is paramagnetic for the 4-leg ladder. 

5) For n p > n (but with n / 1/2 and n not too near 
0), the groundstate of the 2-leg ladder is a paramag- 
netic Luttinger liquid (LL) with charge and spin gaps 
that vanish in the thermodynamic limit to within our 
numerical uncertainty. The ground-states of the 3-leg 
and 4-leg ladders of all lengths we have studied are also 
spin-unpolarized in this range of n, from which we infer 
that the same is true of the 2D system. 

6) For n — 1/2 the 2-leg ladder exhibits a weakly 
dimerized insulating phase with a very small spin-gap. 
However, we have no evidence that this insulating phase 
persists for wider ladders. 

DMRG applied to Hubbard ladders: The Hubbard 
model is defined, as usual, by 
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where cj CT creates an electron with spin polarization a 
on site j and (ij) signifies pairs of nearest-neighbor sites. 
In the limit U — > oo, the Hamiltonian is parameter- free 
in the sense that the second term in H is replaced by 
the non-holonomic constraint of no double-occupancy: 



0, or 1, and one can chose units of energy 



such that t = 1. 
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FIG. 1: The phase diagrams of the infinite U Hubbard model 
on the 2-, 4-, and 6-leg ladders, and the inferred phase dia- 
gram in 2D. 



The DMRG calculations were carried out keeping up 
to 4000^18000 states. All ladders are taken to have open 
boundary conditions in both directions. When we com- 
pute the expectation value of various densities, it is some- 
times useful to break spin rotational symmetry by apply- 
ing a Zeeman field of magnitude h = 1 in the z direction 
on a single boundary site, taken to be the site at the 
lower left-hand end of the ladder. 

To characterize the excitation spectrum of the system, 
we define the charge, spin, and single-particle gaps, A c , 
A s , and Ai p , as follows: 

A c ee [E(N el + 2) + E{N el - 2) - 2E(N el )}/2 (2) 
A s ee [E(S = 1; N e i) - E(S = 0; N e i)} 
A lp ee [E(N el + 1) + E(N et - 1) - 2E(N el )] 

where N e i is the total number of electrons (which is al- 
ways taken, for present purposes, to be even), and E(N e i) 
and E(S,N e i) are, respectively, the ground-state energy 
and the ground-state energy in a given spin-sector. This 
definition of the spin-gap is only useful under circum- 
stances in which the groundstate has 5 = 0. Where pos- 
sible, we have extrapolated values of the gaps to the ther- 
modynamic limit by fitting the data from finite length 
ladders to a quadratic form, A(N) = A + AN^ 1 +BN~ 2 , 
where N is the length of the ladders; unless otherwise 
stated, all gaps listed below have been extrapolated to 
the thermodynamic limit in this way. Details of the var- 
ious extrapolations will be presented at a future date. 

Clearly, in a Fermi liquid (FL), all three gaps vanish 
in the thermodynamic limit. In ID, the FL is unstable 
in the presence of any interactions, but there exists a 
gapless LL which is a direct descendant of the FL. 

Results for the 2-leg ladder: We have computed 
the groundstate properties of the 2-leg ladder as a func- 
tion of n for system sizes 2 x N with N = 20, 30, 40, and 
50. To identify the ferromagnetic portions of the phase 
diagram, we have computed the ground-state magnetiza- 
tion density M — S/S max , where S is the total spin of 
the groundstate, and S max — Nn is the maximum possi- 
ble value of S in a fully spin-polarized state. The results 
are shown in Fig. 3 (a), where different colors denote the 
different system sizes. Since the four curves are nearly 
identical, the extrapolation to the thermodynamic limit 
is trivial. Specifically: 

1) The fully polarized ground state terminates at n = 
n F = 4/5, independent of N [10]. (See Table I.) This 
value of np is not locked by any obvious commensurabil- 
ity effect that we have detected [11]. For instance, if we 
modify the Hamiltonian by making the hopping matrix 
elements on the rungs, t' = 0.5i, where t is the hopping 
matrix element on the legs of the ladder, we find that 
n F = 0.85. 

2) The ground state at n — 3/4 is an insulating para- 
magnetic state with a charge gap, A c = 0.24 ± 0.02t, a 
single-particle gap, A± p = 0.245 ± 0.02i, but a vanishing 
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spin-gap (A s < 3 x 10~ 4 i which is zero within our uncer- 
tainty) . The character of this state on a central segment 
of the ladder is shown in Fig. 2. The thickness of the lines 
on the bonds between sites represents the magnitude of 
the expectation value of the bond-density, 
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(3) 



and the length of the arrow on the site represents the 
expectation value of the spin on that site, 
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Spin-rotation invariance has been explicitly broken in 
this calculation by the application of a Zeeman field on 
the first site on the lower leg of the ladder; translation 
symmetry is broken by the ends of the ladder. There is a 
period-2 bond-density wave, with a magnitude that does 
not decrease with distance from the end of the ladder, 
nor does it depend significantly on the length of the lad- 
der; this signifies a discrete, broken translation symme- 
try. There is also a period-4 ordering tendency of the spin 
density, but with an intensity which decays slowly with 
distance from the end at which the Zeeman field is ap- 
plied. This, and the absence of a spin-gap, signifies that 
there is quasi-long range antiferromagnetic order. What 
is not shown in the figure is that the expectation value 



of the electron density, rij = X) c 



form, 



, is nearly uni- 



n for all j except in the immediate vicinity of 



the ends of the ladder. 

A way to visualize this state is that the ladder sponta- 
neously forms a plaquette state consisting of a checker- 




FIG. 2: Ground-state correlations at n = 3/4 for a central 
portion of (a) the 2 x 20 ladder, (b) the 4 x 20 ladder, (c) 
the 6x8 ladder. The thickness of the lines is proportional 
to the third power of the magnitude of Bij and the length 
of the arrows to the magnitude of Sj , where the axis of spin- 
quantization is set by a Zeeman field of strength h = 1 applied 
to the lower left-hand site of each ladder. The numbers index 
the position along the ladder. For the 2 x 20 ladder, the 
values of Bij in the figure range from Bij — 0.30 (lightest 
line) to Bij = 0.45 (darkest line), while the magnitude of Sj 
ranges from -0.19 to 0.21. For the 4 x 20 ladder Bij ranges 
from Bij = 0.30 to 0.41, while Sj ranges from -0.12 to 0.12. 



board array of weakly coupled plaquettes, each with three 
electrons in a state of total spin 3/2, and a weak, anti- 
ferromangetic exchange coupling between plaquettes. As 
one would expect, the character of this state is not sensi- 
tive to small changes in the Hamiltonian - for instance, 
if we set the rung hopping t' = 0.5< or 2t, we find no 
qualitative change in the insulating phase at n = 3/4. 

3) When n is in the range 4/5 > n > 3/4 the ground 
state appears to be a 2-phase mixture of the fully po- 
larized state with n = 4/5 and the checkerboard state 
with n — 3/4. The evidence for this is as follows: a) As 
seen in Fig. 3 (a), the ground state magnetism is (within 
expected finite size corrections) a linear function of n in 
this range, b) The groundstate energy (not shown) is, to 
similar accuracy, a linear function of n in this range, with 
a continuous first derivative at n = uf; this is precisely 
the behavior expected from a Maxwell construction for a 
two-phase region, c) As a final test, we have applied an 
on-site potential of magnitude /i = 0.3i to the left-most 
2x20 lattice sites of a 2x50 ladder with mean value of 
n = 0.77; as can be seen in Fig. 3 (b), the resulting den- 
sity profile consists of a region with density nj w 4/5 on 
the left portion, and rij ~ 3/4 on the right portion of the 
ladder, with a sharp domain wall separating them. 

4) Again from Fig. 3 (a), it is apparent that for 3/4 > 
n > 3/5, there is a regime in which the ground-state is 
partially spin-polarized, with maximal spin-polarization 
being attained at n = 2/3 where M « 0.5. The behavior 
of the 2-leg ladder in this regime is interesting in its own 
right, but, in contrast to the situation in other ranges 
of n, similar behavior is not seen in wider ladders, as we 
shall see below. We have thus not yet explored the nature 
of the phases that occur in this range of n in any detail. 

5) For 3/5 > n, the groundstate has M = 0. Finite size 
scaling leads to the speculation that, for n ^ 1/2, this is 
a Luttinger liquid phase with A c , A s and A\ p all tending 
to in the thermodynamic limit to within our numerical 




FIG. 3: (a) Magnetization, normalized to its maximum pos- 
sible value, of the 2-leg ladder as a function of n.The black, 
green, blue, and red curves correspond to N=50, 40, 30, and 
For the 6 x 8 ladder Bij ranges from Bij = 0.31 to 0.41, while 20. (b) Expectation value of the site density, rij, in a 2 x 50 



Sj ranges from -0.18 to 0.16. 



ladder with average density n = 0.77 and a chemical potential 
of magnitude [i = 0.3t applied to the left-most 2 x 20 sites. 
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lattice size 


riF 


eo 


lattice size 


riF 


eo 


2x20 


0.800 


-0.6183 


3x20 


0.867 


-0.4442 


2x30 


0.800 


-0.6216 


3x25 


0.867 


-0.4455 


2x40 


0.800 


-0.6233 


3x30 


0.867 


-0.4463 


2x50 


0.800 


-0.6243 


3x35 


0.867 


-0.4469 


4x10 


0.800 


-0.6526 


5x25 


0.827 - 0.833 




4x15 


0.817 


-0.6089 


5x30 


0.820 


-0.6558 


4x20 


0.800 


-0.6613 








6x10 


0.817 


-0.6175 


6x20 


0.808 


-0.6515 


6x15 


0.800 


-0.6733 









TABLE I: The minimum electron density, uf, for which the 
groundstate is fully polarized, M = 1, for various size ladders. 
Here eo denotes the groundstate energy per site when n — uf- 

uncertainty of ~ 0.02i. In the near future we intend to 
perform more extensive studies of this regime to more 
fully explore the nature of the long-range correlations. 

6) For n = 1/2 there is a clearly identifiable charge gap, 
A c ~ O.lt. One can think of this as arising from a state 
in which there is a single electron localized on each rung 
of the ladder, so that the spin-degrees of freedom form 
an effective spin- 1/2 chain, and thus can be expected to 
exhibit one of two possible phases - a gapless phase with 
power law antifcrromagnetic correlations or a dimerized 
phase with a spin-gap. As we shall show in a forthcom- 
ing study, the 2-leg ladder with n = 1/2 exhibits weak 
dimcrization and fairly long-range antifcrromagnetic cor- 
relations, implying that it is in the dimerized phase close 
to the quantum critical point. Further evidence of this is 
provided by studying the changes produced by changing 
the ratio of the hopping matrix elements along the rungs 
and legs of the ladder, t'/t. We find that the dimeriza- 
tion is enhanced and the spin-correlations rendered short- 
range for t'/t = 0.5, while the spin-correlations strength- 
ened for t'/t = 2. 

Results for wider ladders: Wider ladders provide 
necessary clues concerning the evolution of the phase di- 
agram as the 2D limit is approached. However, with in- 
creasing width, it becomes more difficult to obtain fully 
converged results from DMRG; we are thus restricted to 
increasingly shorter ladders as the width increases. 

To study the evolution of the HMF phase, we have 
computed np, the largest value of n for which the ground- 
state is fully spin polarized, for a variety of ladders, with 
the results presented in Table I. An even - odd effect is 
apparent. For the 3 x N ladders np = 0.87, independent 
of N, while for 4 xN (as for 2 xN ) n F = 0.8. For 5- 
and 6-leg ladders, we were restricted to relatively small 
N, but the trend continues, with np slightly greater than 
0.8 for the 5-leg and approximately equal to 0.8 for the 6 
leg. Extrapolating either the even or the odd leg ladder 
results to the 2D limit results in an estimate of n F « 0.8 
as the ladder width tends to infinity. 



We exhibit in Fig. 2(b) and 2(c), respectively, the 
ground-state correlations at n = 3/4 of the longest acces- 
sible 4- and 6-leg ladders. The analogous commensurate 
checkerboard state we found in the 2-leg ladder is mani- 
fest in these correlations, as well. Indeed, the magnitude 
of the broken symmetry does not seem to decrease much 
with increasing ladder width; the strong (weak) bonds 
along legs have typical strengths 0.44 (0.31), 0.38 — 0.41 
(0.30 - 0.31), and 0.36 - 0.41 (0.31 - 0.33) for the 2- 
, 4-, and 6-leg ladder, respectively. The corresponding 
values of the bonds on the strong (weak) rungs of the 4- 
and 6-leg ladder are 0.38 (0.32 - 0.33) and 0.37 - 0.39 
(0.32 — 0.33), respectively. 

For the 4 x 20 ladder with n = 3/4, A c = 0.23±0.006t 
and A s = 0.008 ± 0.002i. For comparison, on the 
2x20 ladder, the charge gap is of comparable magnitude, 
A c = 0.286 ±0.006t, but the spin-gap is more than a fac- 
tor of 20 smaller, A s < 0.0003t. The robustness of the 
charge gap corroborates the existence of a commensurate 
insulating checkerboard state. To the extent that we can 
think of the low energy spin degrees of freedom as cor- 
responding to a spin 3/2 antiferromagnet defined on the 
checkerboard lattice, all (4p + 2)-leg and 4p-leg ladders 
correspond, respectively, to (2p + l)-leg and 2p-leg spin 
3/2 Heisenberg ladders, and accordingly are expected to 
be gapless, or to exhibit a spin-gap, albeit one which 
falls exponentially with increasing n [6] . Together, these 
observations strongly support the conclusion that the 
checkerboard antiferromagnet is the groundstate phase 
in the 2D limit. 

Finally we have found that the groundstates of all the 
4-leg ladders we have investigated have total spin < 2 for 
3/4 > n > 3/5 and total spin < 1 for 3/5 > n respec- 
tively, corresponding to M « 0. We hope to study the 
apparent Luttinger liquid phase of the 2-leg ladder and 
its evolution with ladder width toward the 2D limit us- 
ing a combination of DMRG and variational treatments, 
in order to determine whether there is a well defined 
strong coupling limit [8] of the superconducting state that 
is known to arise at weak coupling [7] . 
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